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Year 2 Progress Report: 

The following efforts were initiated during the last year: 

1. Year-round monthly soil C0 2 flux measurements have been initiated in primary and 
secondary forests and in managed and degraded pastures. In the first year of the project, 
we had only measured in primary forests and in degraded pastures, and measurements 
had only been made four times per year. As described in the attached abstract 
submitted to the Ecological Society of America, this new effort of the second year shows 
the importance of seasonal variation in the managed pastures, where fluxes were very low 
at the end of the dry season (0.09 g m 2 hr' 1 ) and very high in the beginning of the wet 
season (0.21 g m' 2 hr' 1 ). This result differs from our previous work in degraded pastures 
that had relatively low soil C0 2 fluxes in both seasons. As a result of this work, we will 
be developing a new below-ground C budget for the managed pastures, which are 
becoming an increasingly important land-use in this region of the Amazon. 

2. The tedious but essential task of root sorting and weighing has progressed. We now 
have a complete set of samples analyzed for all four ecosystem types at Paragominas, 
and these new data have been used in our most recent model simulations. We have also 
completed analyses of the coarse root fraction at Manaus and Santana; the fine root 
fraction analyses are in progress. 

3. Regional modeling of soil water dynamics and minimum rooting depth has produced a 
surprising result. A Brazilian intern working on this project, Gustavo Negreiros, has 
digitized the RADAMBRASIL soils database and has acquired a 20-year record of 
precipitation for the region. He has developed a "tipping-bucket" hydrologic model that 
is run within a GIS framework. As shown in the attached color figure, the eastern half of 
the Amazon Basin depends on soil water stored >2m depth in order to maintain a 
transpiring (3.5mm/day) evergreen canopy during years of average precipitation. More 
surprising is that, during El Nino years, the greatest requirement for deep soil water 
uptake is centered in northeastern Amazonia, from Santarem to Paragominas in Para 
state, and not at the savanna/forest transition further to the south. Apparently, 
precipitation is most strongly affected in northern Para during El Nino years. A 
simulation of 20% rainfall reduction for the entire basin has the greatest effect at the 
savanna/forest boundary. Hence the type of climate perturbation determines which 
regions are most strongly affected. The results also show that deep rooting must be a 
general phenomenon for much of the basin. 

4. Prototype lysimeters have been designed and installed in soil pits to begin assessing the 
importance of DOC as a source of organic matter in deep soils. Zero-tension lysimeters 
constructed of aluminum have been installed below the litter layer, and they appear to be 
working. Tension lysimeters have also been installed, but we have had problems 
obtaining a sample from the very clayey soils. This will be a focus of the coming year’s 
research. The lysimeter samples will be analyzed for DOC concentration and 14 C content. 
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5. Publications: This has been a fruitful year for our project. The core of our work was 
published in Nature in December, 1994. Three others papers have been published or 
accepted for publication (see list below). The Nature paper and page proofs from Tellus 
are appended. 

Future research: 

The research tasks to be included in the coming year include: 

1. measuring dissolved organic carbon (DOC) throughout the soil profile and estimating 
rates of C input into deep soils via this pathway. 

2. continue year-round monthly soil C0 2 flux measurements in four ecosystems: primary 
forest, secondary forest, managed pasture, and degraded pasture. 

3. revising the model of soil C dynamics for managed pastures. 

4. year-round monitoring of soil moisture dynamics and C0 2 fluxes at a more seasonally 
dry site (Santana de Araguaria) and a less seasonally dry site (Manaus). 


Publications: 

Nepstad, D.C., C.R. de Carvalho, E.A. Davidson, P.H. Jipp, P.A. Lefebvre, G.H. 

Negrieros, E.D. da Silva, T.A. Stone, S.E. Trumbore, and S. Vieira. 1994. The role 
of deep roots in the hydrological and carbon cycles of Amazonian forests and 
pastures. Nature 372:666-669. 

Davidson, E.A., and S.E. Trumbore. Gas diffusivity and production of C0 2 in deep soils 
of the eastern Amazon. Tellus, in press 

Nepstad, D.C., P. Jipp, P. Moutinho, G. Negreiros, and S. Vieira. 1995. Forest recovery 
following pasture abandonment in Amazonia: canopy seasonality, fires resistance, 
and ants. pp. 333-341 In D.J. Rapport, C.L. Gaudet, and P. Calow (Eds.) 
Evaluating and Monitoring the Health of Large-Scale Ecosystems, NATO ASI 
Series, Vol I 28, Springer- Verlag. 

Trumbore, S.E., E.A. Davidson, P.B. de Camargo, D.C. Nepstad, and L.A. Martinelli. 
Below-ground cycling of carbon in forests and pastures of the eastern Amazon. 
Global Biogeochemical Cycles. 

Manuscripts in preparation: 

Camargo, Martinelli, Trumbore, and Nepstad. Soil carbon inventories in mature and 
regrowing forests and pastures of eastern Amazonia. Forest Ecology and Management. 

Nepstad, Vieira, Jipp, de Carvalho. Canopy dynamics, root distribution and soil water 
uptake in deeply-rooting forest and pasture, eastern Amazonia. Oecologia. 

Carvalheiro, K. and D. Nepstad. Deep-soil heterogeneity and fine root distribution in 
Amazonian forest, pasture and secondary forest. J. Tropical Ecology. 
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Nepstad, D. and H. Cattanio. Fine root production to 6 m depth in Amazonian forest, 
pasture and secondary forest. Plant and Soil. 

Nepstad, D., L. Martinelli, P. Camargo, S. Trumbore, E. Davidson, C. Uhl. Complete 
carbon stocks in forest, pasture and secondary forest, eastern Amazonia. 

Negreiros, G., D. Nepstad, E. Davidson, P. Lefebvre, T. Stone. Canopy seasonality, 
available soil water, and the predicted distribution of deep-rooted, evergreen forests in 
Brazilian Amazonia. Global Biogeochemical Cycles 

Stone, T.A., and P. Lefebvre. Analyzing selective logging in eastern Brazilian Amazonia 
over time with Landsat TM data. 

Papers presented (or submitted for presentation) at professional meetings: 

Ecological Society of America, 1994 Annual Meeting, Knoxville TN: 

Davidson et al. Soil carbon dynamics in forests and pastures of the eastern 
Amazon. 

Nepstad et al. Canopy water relations and leaf phenology in deep-rooted forest 
and pasture, eastern Amazonia. 

Ecological Society of America, 1995 Annual Meeting, Snowbird, UT: 

Verchot, L.V., E.A. Davidson, and J.H. Cattanio. Seasonality of trace gas 
emissions from soils in the eastern Amazon. 

Negreiro, G.H., D.C. Nepstad, and E.A. Davidson. Mapping the minimum rooting 
depth of Amazonian forests. 

Nepstad, D.C., P.R. Moutinho, T. Restom, and V. Mendes. Deep-soil biology of 
an Amazonian forest: mycorrhizae, soil heterogeneity, and cutter ants. 

First GCTE Science Conference, May 1994, Woods Hole, MA: 

Davidson et al. Soil carbon dynamics in forests and pastures of the eastern 
Amazon. 

Nepstad et al. Soil water extraction, canopy seasonality, and drought-stress in deep 
rooted primary forest, pasture and secondary forest, east Amazonia. 

International Symposium on Resource and Environmental Monitoring Sept., 1994, Rio de 

Janeiro: 

Negreiros et al. The distribution of evergreen forests in seasonally-dry Amazonia. 

Lefebvre P.A. and T.A. Stone. Monitoring selected logging in eastern Brazilian 
Amazonia using multi-temporal Landsat Thematic Mapper imagery. 

Soil Science Society of America, 1994 Annual Meeting, Seattle WA: 

Trumbore, S.E. and E.A. Davidson. Gaseous diffusivity and 

production of C0 2 in deep forest soils of the Eastern Amazon 
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MINIMUM ROOTING DEPTH 
WITH AVERAGE RAINFALL 


Preliminary results of Negreiros et al.: 


Minimum rooting depth of 
evergreen forest in Amazonia, 
calculated within a geographical 
information system based on a 
regional map of soil water holding 
capacity (calculated from 1140 soil 
profiles, soil water retention 
equations from the literature, 
rainfall data from 600 weather 
stations, and an estimated 
evapotranspiration of 3.5 mm/day 
(Nepstad et al., 1994). These 
versions of the map consider 
three scenarios: average monthly 
rainfall, the rainfall pattern of an 
ENSO year (1983), and a 20% 
reduction in rainfall. The study 
sites for the research of this 
project are indicated in the 
second map. 



MINIMUM ROOTING DEPTH 
DURING ENSO (1983) 


Trombetas 



MINIMUM ROOTING DEPTH 
WITH 20% RAINFALL REDUCTION 



Min. roofing depth 
0 - 2 m 

i 2 - 4 m 

iH > 4 m 


Savanna and Deciduous Forests 
Cleared for agriculture 
Seasonally Flooded and other 
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LETTERS TO NATURE 


The role of deep roots in the 
hydrological and carbon cycles 
of Amazonian forests 
and pastures 

Daniel C. Nepstad*, Claudio R. de Carvalho' , 

Eric A. Davidson*, Peter H. Jlpp*j, 

Paul A. Lefebvre*, Gustavo H. Negrelros*, 

Elson D. da Sllvat, Thomas A. Stone*, 

Susan E. Trumbore§ & Simone Vieira* 

USA°° dS H °' e Research Center ' Woods Hole, Massachusetts 02543, 

t EMBRAPA-CPATU, CP 48, Belem, Para 66.001, Brazil 
§ Department of Earth System Sciences, University of California 
Irvine, California 92717, USA 


Deforestation and logging transform more forest in eastern and 
southern Amazonia than in any other region of the world 1 "’. This 
forest alteration affects regional hydrology 4 " 11 and the global car- 
bon cycle " , but current analyses of these effects neglect an 
important deep-soil link between the water and carbon cycles. 
Using rainfall data, satellite imagery and field studies, we estimate 
here that half of the closed forests of Brazilian Amazonia depend 
on deep root systems to maintain green canopies during the dry 
season. Evergreen forests in northeastern Para state maintain 
evapotranspiration during five-month dry periods by absorbing 
water from the soil to depths of more than 8 m. In contrast, 
although the degraded pastures of this region also contain deep- 
rooted woody plants, most pasture plants substantial^ reduce their 
leaf canopy in response to seasonal drought, thus reducing dry- 
season evapotranspiration and increasing potential subsurface run- 
off relative to the forests they replace. Deep roots that extract 
water also provide carbon to the soil. The forest soil below I m 
depth contains more carbon than does above-ground biomass, and 
as much as 15% of this deep-soil carbon turns over on annual 
or decadal timescales. Thus, forest alteration that affects depth 
distributions of carbon inputs from roots may also affect net carbon 
storage in the soil. 

Pastures are the most common type of vegetation on de- 
forested land in Amazonia. They vary greatly in the ratio of 
grass to woody-plant cover, including managed pastures (from 
which woody vegetation is removed by heavy machinery) and 
the more common shrub- and tree-dominated “degraded" 
pastures . We studied deep roots in a mature, evergreen forest 
and in adjacent pastures near the town of Paragominas, in the 
Brazilian state of Para. The deeply weathered clay soils at the 

L^ rt , 7 a r^i ,ePartment ° f En '' ironmemal Studies. Duke University, Durham, North 
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FIG. 1 Major forest types and dry-season rainfall 
of Brazilian Amazonia. Savannas and deciduo 1 ^ 
forests (14%) were separated from evergreen 
forests (75%) based on seasonal patterns of 
canopy greenness (as seen from satellites) and a 
vegetation map. Evergreen forests include areas 
that did not display a seasonal depression of the 
normalized difference vegetation index 
(NDV!) ?8?9 during the dry seasons of 1986-88. 
To minimize the effects of clouds and smoke, 
only maximum monthly values of NDVI were 
used. Areas cleared for agriculture and season- 
ally flooded land (that is, pixels with the spectral 
signature of water) were adapted from Stone et 
a/ 29 . Dry-season rainfall refers to the driest 
three-month period of the year averaged for each 
rainfall record 30 Of the 212 weather stations 
used, 75% had more than eight years of rainfall 
data. Interpolation of rainfall data was conducted 
using a Geographic Information System software 
package (IDRISI; Clark Univ., Worcester, Massa- 
chusetts). Annual rainfall is 1,500 mm at San- 
tana de Araguaia and daily rainfall is <0.5 mm 
during the driest three months. Annual rainfall at 
Manaus and Trombetas is 2,300 mm and daily 
rainfall is >1.5 mm during the driest three 
months. Evergreen forests in areas with 
<1 5 mm d 1 during the driest three months of 
the year cover 36% of the region, and must 
depend on soil water extraction below 1 m depth, 
as described for the Paragominas study site 
(2 59' S, 47 31’ W). 



site are common in Amazonia’". Precipitation is highly seasonal 
with an average or 1.750 mm annually and <250 mm trom July 
to November. This seasonality is typical of the eastern and sou- 
thern portions of Amazonia where most deforestation is occur- 

rll1 Root^ were most abundant near the soil surface, as expected 
but were found in deep shafts to depths of ~8 m in the managed 
pasture. -12 m in the degraded pasture, and -18 m in the foresU 
Fine-root biomass was lower in the top 10 cm of the degraded 
pasture soil than in the forest, but no differences occurred at 
greater depths. In the managed pasture, fine-root biomass was 
three times greater than in the forest near the soil surface and 
15 100 times lower at depth (Fig. 2). 

We assessed the hydrological role of deeply penetrating ro 
using a soil water-balance approach. During the severe 5.5- 
month dry season of 1992, when total rainfall was only 95 mm, 
plant-available soil water (PAW) from 2-8 m depth dee med by 
580 mm in the forest and 310 mm in the degraded pasture. In 
the upper 2 m of soil, the decline in PAW for _ this period 1 was 
130 mm and 100 mm, respectively (Fig. 3a. b). Hence PAW 
stored below 2 m in the soil provided >75% of the water 


FIG 2 Vertical profile of live fine-root biomass (diameter <1 mm) in 
adjacent mature forest and man-made pastures near Paragominas, 
eSem Amazonia. Points are means, and bars are one standard error 
of the mean from 1.5-kg soil samples taken at each ® ar ^P lin ® d |^ h h 
36 auger borings in the forest, 17 in the degraded pastures, and 6 in the 
managed pasture. Roots were cleaned by flotation sieving and sorted 
manually at lOx magnification. The forest was cleared at the P aatL | re 
sites in 1969, planted with the grasses Panicum maximum and later 
Brachiaria humidicola. and was heavily grazed '^te mitten y 
present. Woody shrubs and treelets now comprise 50% of live leaf area 
in the degraded pastures and they support little grazing. The ™ ana S e d 
pasture was disk-harrowed, fertilized and replanted with Brachiaria bri- 
zantha in 1988, and is virtually free of woody vegetation. 


traded from the soil in both ecosystems, and 100 mm more 
iter was extracted from the forest soil. This is the case despite 
e fact that the deep forest soil began the 1992 dry season 
ithout being fully recharged (Fig. 3 b) because ol below-average 
linfall in the wet season. More rainfall is needed to recharge 
le soil of the forest because it extracts more soil water during 
,e dry season and because its canopy intercepts a larger portion 

f incoming rain. , r 

Differences in water extraction between the forest and 
esraded pasture result in contrasting patterns of evapotranspir- 
tion Evapotranspiration can be estimated for the dry season 
•om the sum of rainfall and PAW depletion (surface runoff is 
isignificant in both ecosystems because of high soil infiltration 
ates 19 ). During the dry season of 1992. average daily rain ( a1 ’ 
/as 0.6 mm, and evapotranspiration in the forest and degraded 
iasture was 3.6 and 3.0 mm, respectively. Current vegetation - 
tmosphere models developed to predict the effect on climate of 
Amazonian deforestation assume rooting depths of only 1.3- 
(.0 m in forests and 0.6 m in pasture 8 and, therefore, would 

inderestimate evapotranspiration by >60% in these Para- 
i inm rit-ii enacfin Important 


Fine-root biomass (mg cm" 3 ) 
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inerr 0 ’- 1 . But the soil C inventory below 1 in depth at Para- 
gominas exceeds soil C above 1 m (Fig. 4 a), as well as exceeding 
above-ground biomass ( 18 kg C in 2 ). Contrary to the assump- 
tKm of inertness, our isotopic data indicate the presence of a 
significant fraction of modern C in deep soils. First, the A l4 CO, 
(sec Fig. 4 legend for definition of this quantity) below 1 m 
depth remains at or above modern atmospheric levels (Fig. 4b ), 
indicating that the high concentrations of CO : observed at depth 
(Fig. 4c) result from root respiration and from microbial decay 
of carbon fixed by plants within the past 30 years 22 . Second, the 
A C of soil organic matter declines with depth, as expected, 
because the proportion of very old C increases with depth (Fig. 
4^/). but 10 15% of the soil C at 8 m depth may be modern (if 
the old C is essentially radiocarbon dead at a A 14 C value of 
-l,000% o , the remaining modern fraction at +143% 0 to +200%o 
would be -13% of the total C present at 8 m depth 
[( -850+ 1 ,000)/ 1 ,000]/[ 1 +( 1 43/1 ,000)] =0. 1 3). This result sug- 
gests that up to 3 kg C m 2 occurring below 1 m depth cycles on 
annual to decadal timescales, and that this C would be subject to 
change if root distributions were changed by land-use practices. 

In the degraded pastures, woody plants have deep-root sys- 
tems (Fig. 2) that may maintain the deep-soil C pool if fine- 
root turnover rates arc similar to those in the forest. We have 
concentrated our studies of soil C dynamics on forests and 
degraded pastures, as these ecosystems are clearly the most 
abundant at present in eastern Amazonia. However, disk-har- 
rowing and herd rotation can effectively exclude woody vegeta- 
tion in more intensively managed pastures, and these practices 
are becoming more common 1617 . Based on root distribution 
(Fig. 2), we would predict that these managed pastures would 
lose deep-soil C and possibly gain surface-soil C (a hypothesis 
which we now intend to test). Differences in pasture management 
that affect distributions of C inputs throughout the soil profile 
may help explain why some studies have showm increases while 
others have shown decreases in C inventories of the surface soil 
following conversion of forest to pasture 23 25 . 

How common is deep rooting in the Amazon Basin? Roots 
have been inferred to 5 m depth in a forest of Surinam 27 . Our 
shaft excavations to 8 m depth at forest sites that are less sea- 
sonal (Trombetas, Manaus) and more seasonal (Santana de Ara- 
guaia) than Paragominas (Fig. 1) revealed root distributions 
very similar to that of the Paragominas forest. These initial field 
data indicate that deep rooting is common in Amazonia. 

We estimated the geographical distribution of deeply rooting 
forests in Brazilian Amazonia by overlaying monthly estimates 
of canopy greenness from AVHRR satellite imagery with a Geo- 
graphic Information Systems database of rainfall from 212 
weather stations. Where forests are evergreen but seasonal 
drought is significant ( <1.5 mm d 1 during the driest three 
months) we deduce that the forest must rely on water uptake 
from deep soil. The forests meeting these criteria cover an area 
of -1.8x JO 6 km 2 , which is most of the eastern and southern 
half of the Amazonian closed-canopy forest (Fig. 1 ). Hence deep 
roots play an important role in maintaining dry-season canopy 
greenness and evapotranspiration in the regions where human 
activity is concentrated. Deep roots are not limited to seasonally 
dry regions and may also play a role in nutrient uptake. 

Deep roots help explain why Amazonian evergreen forests 
extend well into a region characterized by a long dry season. 
Understanding the effect of human land-use practices on 
regional budgets of water and carbon will require knowledge of 
the basic processes involving deep roots and deep soils. □ 


8. Nobre. C., Sellers. P. & Shukla, J. J. Clim. 4 , 957-988 (1991). 

9. Shukla. J.. Nobre, C. A. & Sellers. P. Science 247 , 1322-1325 (1990). 

10. Lean, J. & Warrilow, D. A. Nature 342 , 411-413 (1989) 
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hydrological processes, involving seasonal variation of deep 
water storage as high as 400 mm (Fig, 3b ), are left out of models 
that ignore deep roots. 

Evapotranspiration may have been reduced in the degraded 
pasture by a loss of green leaf area. The leaf area of degraded 
pasture vegetation declined 68% whereas average leaf area o 
forest plants declined only 16% during the 1992 dry season (Fig. 
3c’ ; the managed pasture lost 100% of its leaf area during this 

same period). , 

Less depletion of PAW in the degraded pasture signifies that 
this ecosystem can store less rainfall than the forest, and may 
therefore produce far more seepage to the groundwater aquifer 
or subsurface runoff to streams during the wet season. At the 
end of the 1992 dry season, the forest could store an additional 


FIG. 3 Seasonal trends in soil water and leaf area in adjacent mature 
forest and degraded pasture, eastern Amazonia. Points are means, and 
bars are one standard error, a, Plant-available soil water (PAW) from 
0-2 m depth, b, PAW from 2-8 m depth, c, Percentage of maximum 
leaf area of ten common species in the forest and degraded pasture, 
based on monthly observations of tagged branches (n = 20 branches 
per species). In degraded pasture, values were weighted by the relative 
abundance of grasses and non-grasses, d, Daily rainfall. PAW was cal- 
culated as the amount of soil water held at tensions between 
-0.01 MPa (field capacity) and -1.5 MPa, determined from soil mots- 
ture retention curves for intact soil samples taken at various depths. 
Biweekly measurements of volumetric soil water content were made 
using time domain refiectometry, with sensors (6-8 per depth) installed 
at the surface and 1.5 m horizontally into the shaft walls at 1-m depth 
intervals 31 ' 32 . 


770 mm of water in the upper 8 m of soil, compared to 400 mm 
in the pasture (Fig. 3a, b). 

We assessed the role of root distributions in the soil carbon 
cycle by assessing soil C stocks and by analysing the isotopic 
composition of soil C0 2 and soil organic matter. The degraded 
pasture soil has lost 1.6 kg C m~ 2 in the top 1 m of soil relative 
to the forest (Fig. 4 a); this loss is coincident with lower fine- 
root mass (and presumably lower C inputs) in the surface soil 
of the degraded pasture (Fig. 2). Changes in soil C stocks at 
depth are more difficult to determine reliably because small 
differences in C concentrations that are near detection limits 
must be extrapolated over large soil volumes. Soil C inventories 
and models of soil C processes seldom include analyses below 
1 m, and the C present in deep soils is often assumed to be 


Soil C (kg m* 2 ) 




A 14 C0 2 (%*) 



A 14 C in SOM ( %o) 

\ Q00 -800 -600 -400 -200 0 




1992 

ATMOSPHER" 




g Degraded 
pasture 


FIG. 4 Vertical profiles of soil C properties; points are means and bars 
are one standard error of the mean of three soil shafts per ecosystem, 
a Soil C stocks (not including root C). b, A 14 C0, of the soil atmosphere, 
c! Concentration of C0 2 in the soil atmosphere, d, A C of soil organic 
matter (SOM) excluding root C. Soil gases were sampled through stain- 
less-steel tubes (1.8 m long, 0.32 cm diameter) inserted horizontally 
into the sides of the soil shafts. Analysis of radon activities showed that 
the effect of gas exchange through the shaft wail was negligible. C0 2 
concentrations were analysed in the field using a LiCor (Logan, Utah) 
infrared gas analyser. Gas samples were also stored in 0.5-1 stainless- 
steel canisters for later analysis of 14 C0 2 by accelerator mass spectro- 
metry. Reported as A 14 C (the 14 C/ 12 C ratio in parts per thousand cor- 
rected for variations in 13 C/ 12 C and compared to an absolute standard 
of oxalic acid in 1950), zero represents the 1950 atmospheric C0 2 
content, positive values reflect the influence of C derived from above- 
ground thermonuclear weapons testing, and negative values indicate 
that significant portions of the 0 have been isolated from 
atmospheric C0 2 long enough for the C to decay (half-life, 5,730 

years) 22 . 
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ABSTRACT 

The rate of gaseous diffusion in soils affects the exchange of gases between the soil and the 
atmosphere, thereby affecting rates of soil respiration and other soil microbial processes. 
Understanding the causes of spatial and temporal variation in soil diffusivity will help explain 
controls of soil sources and sinks of atmospheric gases. In a study of sources of C0 2 in deep soils 
of forests and pastures of the eastern Amazon, we estimated gaseous difFusivity from bulk density 
and volumetric water content using published equations that assume the soil to be either an 
aggregated or nonaggregated medium. The aggregated model requires differentiation of inter- 
and intra -aggregate pore space; we estimated intra-aggregate pore space from volumetric water 
content at field capacity. Steady state m Rn profiles were predicted from a 1-D model using the 
difiusivities generated by both aggregated and nonaggregated models. Predicted values were 
compared with M2 Rn activities measured to 5 m depth- While the models predict similar radon 
activities below about 1 m, large differences are predicted for the top 1 m of soil. The non- 
aggregated model underestimated diffusivity and overestimated ^Rn activities at 1 m and 
above, which is not surprising given that surface soils are usually well aggregated. Having 
validated the aggregated media model using the 222 Rn profiles, estimates of diffusivity were 
combined with measured profiles of C0 2 concentrations to estimate C0 2 production by 
depth. About 70-80% of the measured C0 2 flux from the soil surface was produced in the top 
1 m of soil (including litter in the forest). The 20-30% produced below 1 m results from root 
respiration and microbial decay of root inputs at depth, indicating that deep soil processes are 
a non-trivial component of carbon cycling in these deep-rooting ecosystems. About 1 % of the 
20 kg C m ~ 2 stock of soil C found between 1 m and 8 m depths turns over annually, indicating 
that land-use changes that affect rooting depth could significantly affect deep soil C stocks over 
decades to centuries. Fully understanding the role of land-use change on the global carbon cycle 
will require consideration of these deep soil processes. 


1. Introduction 

Because soils store two to three times as much 
carbon as the atmosphere, the effect of changing 
climate on soil respiration and other soil micro- 
bial processes is important as a potential feedback 
mechanism of global warming (Billings, 1987; 
Oechel et al., 1993; Woodwcll, 1989). Deforesta- 
tion and other land-use changes also affect the soil 
environment and contribute to transfer of C from 


* Corresponding author. 


soils to the atmosphere (Davidson and Ackerman, 
1993). 

Production of C0 2 by soil microorganisms 
roughly equals terrestrial net primary productivity 
(NPP) on a global basis (give or take a couple 
of gigatons of C per year; Schlesinger, 1991). 
Above-ground processes such as NPP are more 
easily studied than are below ground processes 
such as soil respiration. For example, relatively 
sophisticated physiological models, based on 
first-principal understanding of the biochemical 
mechanisms of plant photosynthesis and respira- 


Tellus 47 B (1995), 0 



002 


E. A. DAVIDSON AND S. E. TRUMBORE 


tion, have been developed in so-called “big-leaf 1 
models at local and global scales (Amthor, 1994; 
Sellers etal., 1992). In contrast, the best soil 
respiration model may be a simple temperature- 
dependent Q l0 function based on empirical fits of 
data (Raich and Schlesinger, 1992). Soil is difficult 
to model because it is a complex medium that con- 
sists of a broad range of types of organo-mineral 
particles and aggregates and that contains 
numerous organisms with differing physiological 
processes. Soil properties vary temporally and 
spatially, both horizontally and vertically. 

Although the complexity of the soil matrix often 
seems daunting, some soil processes are under- 
stood well enough to make progress towards 
characterizing probable responses of soils to 
changes in temperature, precipitation, and land 
use. One such area of research is the influence of 
rates of diffusion of gases on soil microbial activity 
and on fluxes of trace gases such as C0 2 , CH 4 , 
N 2 0, and NO (Billings, 1987; Davidson, 1993; 
Born et al., 1990; Oberbauer et al., 1992; Roulet 
ct al., 1992). The presence of 0 2 is an important 
controller of rates of nitrification, denitrification, 
methane oxidation, methanogenesis, and aerobic 
respiration (Conrad, 1989; Firestone and David- 
son, 1989; Linn and Doran, 1984; Schimel et al., 
1993). Understanding the factors that affect diffu- 
sion of 0 2 into the soil and diffusion of trace gases 
out of the soil is necessary to describe variation in 
soil respiration and other soil microbial processes. 

The diffusion of gases through soil is dependent 
on soil porosity and soil water content. Direct 
measurements of diffusivity of gases in soils are 
difficult and cumbersome (Dorr and Miinnich, 
1990; Rolston etal., 1991). Alternatively, gas dif- 
fusivity can be calculated from models that require 
estimates of soil porosity and soil water content. 
Because it is easier to monitor changes in soil 
water content than to make frequent direct 
measurements of effective diffusivity, use of models 
offers an obvious advantage for studying the 
dynamics of gas diffusivity in soil. 

The models that have been proposed for 
estimating gas diffusivity range in complexity, data 
requirements, and usefulness. The most simple 
models (Millington, 1959; Millington and Quirk, 
1961 ) were intended for nonaggregated media, but 
have been used, nevertheless, for soils (Hendry 
et al, 1993; Wood et al., 1993). A model for aggre- 
gated media (Millington and Shearer, 1971) is 


more complicated and requires estimates of inter- 
aggregate and intra-aggregate porosity and water 
content. DilTusivity can also be estimated in soils if 
the distribution of pore sizes can be estimated 
(Nielson et al., 1984). 

The purpose of this research was to develop and 
test a model of soil gas diffusivity that estimates 
C0 2 production as a function of depth using only 
input data on soil profiles of bulk density, C0 2 
concentrations, ambient water content, and water 
content at field capacity. We are interested in C0 2 
production deep in the soil because the deeply 
weathered Oxisols of the eastern Amazon basin 
contain >20 kg C m“ 2 at depths below 1 m, and 
roots extend to at least 18 m depth (Nepstad et al., 
1994). The turnover time of C in deep soils is not 
well known, and study of C0 2 production is one 
means of addressing the dynamics of this deep soil 
C. We measured profiles of 222 Rn, a radioactive 
noble gas (half-life of 3.8 days) produced in soils 
by decay of 226 Ra, to test models of soil diffusivity. 
We then applied the model of Millington and 
Shearer (1971) for aggregated media to measured 
profiles of C0 2 and water content to show how 
C0 2 production varies with depth between wet 
and dry seasons and among forest and pasture 
sites. 


2. Background theory 

One of the earliest and most commonly used 
models of gaseous diffusion in soil (Penman 1940) 
is a simple linear relation: 

DJD 0 = xa , (1) 

where Z) s is the diffusion coefficient in soil, D a is 
the diffusion coefficient in air (e g., 0.135 cm 2 s -1 
for 222 Rn and 0.162 cm 2 s " 1 for C0 2 at 25°C and 
standard pressure), x is a dimensionless tortuosity 
factor, and a is the air-filled porosity (cm 3 air 
cm -3 bed space). Penman (1940) reported a value 
of 0.66 for t. However, tortuosity can also be 
a function of soil water content (Millington, 
1959), because the pathway for gaseous diffusion 
becomes more tortuous as the soil becomes wetter. 
Hence, t must be derived for each soil water con- 
tent of interest for each soil of interest. 

Millington (1959) proposed the following rela- 
tion for diffusion of gases through non-aggregatcd, 


Tell us 47B (1995), 0 




003 


GAS D1FFUSIVITY AND PRODUCTION OF C0 2 


dry, porous media that is assumed to consist of 
interpenetrating spherical solids separated by 
interpenetrating spherical pores: 

DJD 0 — a 4/3 . (2) 

In dry media, a is equivalent to total porosity. As 
the porosity decreases (due to compaction of a 
given medium or differences among types of media 
compared) the probability of pore spaces being 
continuous declines nonlinearly. 

For wet media, the probability of continuous 
pore spaces declines more rapidly with decreasing 
air-filled porosity than accounted for in eq. (2). 
Millington (1959) added a second factor to the 
equation so that diffusivity is affected not only by 
the amount of air-filled pore space as in eq. (2), but 
also by the square of the fraction of total pore 
space that is air- filled: 

DJD 0 = a 4 '\a/e)\ (3) 


into intra-aggregate pore space (e ri ; cm 3 intra- 
aggregate pore space cm -3 bed space) and inter- 
aggregate pore space (e w ; cm 3 inter-aggregate 
pore space cm " 3 bed space). The equation has two 
terms: the first is for diffusion within the intra- 
aggregate space, and the second is for diffusion 
within the inter-aggregate space. The second term 
is identical to eq. (4), except that the parameters 
refer specifically to total porosity and air-filled 
porosity of the inter-aggregate space only: 


DJVo 





(7) 


where i is the total porosity (cm 3 pore space cm 3 
bed space). 

In a further refinement of this model, the 
exponent of the first term of eq. (3) varies to 
improve the statistical function for the probability 
of continuity of pore space within the medium 
(Collin and Rasmuson, 1988; Millington and 
Quirk, 1961): 

DJD 0 = a 2 '(a/e)\ (4) 

where x is determined from the relation: 

a 1 ' + (1 —a) x = 1. (5) 

The value of x must be solved numerically and 
typically ranges from 0.6 to 0.8 (Collin and 
Rasmuson, 1988); 2x is usually near 4/3. The value 
for x can also be approximated from the poly- 
nomial: 

.v = 0.477a 3 - 0.596a 2 + 0.437a 4- 0.564. (6) 

It is important to note that eqs. (2)— (5 ) are 
intended for application to non-aggregated media. 
Soils are usually aggregated, so it would not be 
surprising if these equations do not accurately 
predict diffusion of gases through soils. Millington 
and Shearer (1971 ) have proposed a related model 
for diffusion through aggregated media, but it 
requires that the soil be conceptually divided 


. where a„ is the air-filled intra-aggregate pore 
space (cm 3 air cm -3 bed space), a er is the air-filled 
! inter-aggregate pore space, S is the space occupied 
by solid (cm 3 solids cm" 3 bed space), and where x, 
y , and z are derived from the following equations 
which are similar to eq. (5): 



# + (l 1 , 
°lr + (1 — a tiY = 1 * 


( 8 ) 

(9) 

( 10 ) 


The length and seemingly complicated nature of 
eqs. ( 7 )— ( 10) may discourage some researchers 
from attempting to use them. More importantly, 
estimates of intra- and inter-aggregate porosity are 
often not available. On the other hand, using 
a mode! designed for aggregated media, which 
most soils clearly are, may provide a significant 
improvement over use of models designed for non- 
aggregated media, even if only crude estimates of 
the division between intra- and inter-aggregate 
porosity are available. 

One approach for estimating intra- and inter- 
aggregate porosities is to assume that all intra- 
aggregate spaces, and no others, are entirely 
water-filled at field capacity and, hence, can be 
estimated from volumetric water content at field 
capacity (0 fc ). Field capacity is defined as “the 
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content of water, on a mass or volume basis, 
remaining in a soil 2 or 3 days after having been 
wetted with water and after free drainage is negli- 
gible” (Soil Science Society of America, 1987). The 
inter-aggregate pore space is assumed to be the 
macropores that are water-filled only when the soil 
is saturated, and they are air-filled at field capacity. 
Hence, in ter- aggregate porosity is estimated 
from the difference between total porosity and 
volumetric water content at field capacity (s — 0 fc ). 
Sollins and Radulovich (1988) observed rapid 
drainage of saturated volcanic soils of Costa Rica 
until the matric potential reached —0.05 MPa and 
very slow drainage thereafter. Their soils were 
highly aggregated, and they deduced that the 
difference between water content at saturation 
(when all pore spaces are water-filled) and at field 
capacity (which they defined as —0.03 MPa for 
their soils) provided an estimate of the inter- 
aggregate porosity. Their purpose was to study 
biphasic flow of water rather than gases, but the 
same principals should apply to either fluid. 

3. Methods 

3.1 . Site description 

The study was conducted at the Fazenda Vitoria 
(Victory Ranch), located 6.5 km northwest of the 
town of Paragominas, Par& State, Brazil, in east- 
ern Amazonia (2°59'S, 47 0 31'W). Paragominas 
has been a center of cattle ranging and logging 
since the 1960’s (Nepstad et al., 1991). Average 
annual rainfall is 1750 mm and is highly seasonal; 
<250 mm fall from July to November in an 
average year. The forest maintains an evergreen 
canopy throughout the dry season, apparently 
by extracting water from deep in the soil profile 
(Nepstad et al., 1994). The soils are deeply 
weathered, well-drained, kaolinitic, red-yellow 
Oxisols (Haplustox). This climate is typical and 
these soils are common within the “arc of 
deforestation” ranging from the northeast to the 
southwest of the Brazilian Amazon (Nepstad et al., 
1991 ). The water table is at about 45 m depth. 

3.2. Instrumentation of soil pits 

Soil pits were dug by hand in forest and pasture 
study plots. The pit openings were I m x 2 m and 
they extended to 9 m depth. Three soil pits in 
forests and three in degraded pasture were studied. 


To keep illustrative figures from being cluttered, 
representative data from only one pit per eco- 
system type are shown in figures, but data from all 
pits were used to analyze radon profiles and to 
calculate depth distribution of C0 2 production. 

Pairs of time domain reflectometry (TDR) 
probes were installed horizontally into the pit wall 
at 50, 100, 200, 300, 400, 500, 600, 700, and 800 cm 
depth. Holes (10 cm diameter) were augered 
150 cm into the pit walls, the 30 cm long probes 
were then further inserted into undisturbed soil at 
the end of the auger holes, and the holes were 
backfilled. A pair of probes was also installed verti- 
cally into the soil surface between 0 and 30 cm 
depth at each study site. The study sites have 
been visited weekly for TDR measurements since 
1991. The TDR measurements were converted to 
volumetric water content using calibrations curves 
derived from laboratory analysis of intact soil 
cores. The calibration curves are nearly identical to 
published calibrations of TDR for a variety of 
other soils (Topp et al., 1980; Topp and Davis, 

. 1985). 

Stainless steel tubes (3 mm OD) were also 
installed horizontally into the pit walls at about 15, 
25, 50, 75, 100, 300, 500, and 800 cm depth, and 
the auger holes (10 cm diameter) were backfilled, 
except that a chamber about 10 cm long was left 
open at the far end of the tube buried in the soil. At 
100 cm depth and lower, the tubes were 180 cm 
long. Above 1 m depth the tubes were 90 cm long. 
Perforations were made in the last 10 cm of tubing 
to act as inlet vents. The ends protruding from the 
pit walls were fitted with Swagelok fittings and 
septa to permit sampling of soil gases with a needle 
and syringe. 

3.3. Measurement of CO 2 flux from the soil surface 

Flux measurements were made by circulating air 
between a LiCor infrared gas analyzer and a 
flux chamber consisting of a PVC ring (20 cm 
diameter x 10 cm height) and a vented PVC cover 
(10 cm height). In the forest, the rings were placed 
over the forest floor so as not to disturb the litter 
layer, and wetted soil that is very low in carbon 
content (taken from deep in the soil pits) was 
packed around the outside of each ring to main- 
tain a seal; this deep soil produced a small pulse of 
C0 2 when it was mixed, but C0 2 emissions were 
negligible after one day, making it suitable for 
sealing the outside of the chambers. In the pasture 
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where there is no litter layer, the rings were 
inserted 1 to 2 cm into the mineral soil. Green 
plant tissue was clipped prior to placing the cover 
over each ring. The chamber was vented to the 
atmosphere through a stainless steel tube (10 cm 
long x 3.2 mm OD) to keep pressure within the 
chamber equal to the atmosphere. A battery- 
operated pump maintained a flow of 0.7 1/min to 
and from the chamber and the LiCor. Varying the 
flow rate from 0.4 to 1.2 1/min had no detectable 
effect on measured flux rates. The LiCor readout 
was recorded on a 15-s interval with a datalogger 
and the C0 2 flux was calculated from linear regres- 
sion of the increasing C0 2 concentrations within 
the chamber between 1 and 5 min after placing the 
cover over the ring. 8 flux measurements were 
taken per site and date; coefficients of variation 
were typically 30%. The LiCor was calibrated 
in the field by allowing the pump to draw air 
from a vented tube through which a standard of 
377 ppm C0 2 in air was flowing, while using the 
same lengths of tubing as in flux measurements 
to ensure that the LiCor experienced similar pres- 
sure differentials during flux measurements and 
calibrations. 

3 A. Soil CO 2 concentration depth profiles 

The stainless steel tubes inserted into the pit 
walls were first flushed by removing a 20 ml 
sample with a syringe and discarding that sample. 
Duplicate 3 ml samples were then drawn from 
each tube and the syringes were closed with nylon 
stopcocks and were carried to the soil surface. 
A 1 ml sample was then immediately injected from 
these syringes into a stream of C0 2 -free air passing 
through the LiCor (Fig. 1 ). The C0 2 -free air was 
produced by filtering air drawn from a battery 
operated pump through a soda lime scrubber. As 
the pulse of syringe air drawn from the soil passed 


through the LiCor, a peak was recorded with a 
datalogger using 1 s recording intervals. Injections 
of 0.5 to 3.0fotf**of a standard (2.5 % C0 2 in air) 
were also injected to produce a calibration curve 
: relating peak area to C0 2 molar volume. The 
coefficient of variation for replicate injections of 
standard gases was usually <2%. The calibration 
gas concentration was certified to 2%. 

To check for the effects of removing large 
volumes of gas from the soil on the concentration 
of C0 2 measured, we removed five 20 ml samples 
from a single gas sampling tube and analyzed a 
1 ml subsample from each. After removing a total 
of 100 ml of soil gas from a single tube, we found 
no change in C0 2 concentration. 

5,5. Soil radon activity depth profiles 

222 Rn measurements were made by direct alpha 
counting of air samples. Samples of 50 ml were 
j drawn from soil sampling tubes, dried (by passing 
j over a column of Drierite, which is anhydrous 
j calcium sulfate) and introduced into pre-evacuated, 
100 ml phosphor-coated Lucas counting cells 
(Mathieu et al., 1988; Trumbore et al, 1990). 
Ambient air was used to backfill the cells to atmo- 
spheric pressure. The cells were stored for at least 
3 hours after filling to establish secular equilibrium 
j of 222 Rn with two short-lived, alpha-emitting 
daughters. Scintillation counting of gross alphas 
(using a PylonR counter) was performed within 
. 24 h of sample collection. Radon activities were 
calculated according to the formula: 

C-B 

^“ ,r_ axCEx Kxexp <- ' i( ' n,p ~ lf<:,> ’ 

where A cxu is the radon activity at the time of 
sample extraction from the soil; C is the number of 
counts per minute measured by the counter; B is 


ambient 

air — ► A3 


pump 


soda lime 


scrubber 
(120 cc) 



Fig. 1. Diagram of apparatus for field measurement of COj concentrations of syringe samples of soil gases. 
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Table 1. Soil properties , volumetric soil water content , and calculated effective diffusivities for C0 2 by depth in a forest and a pasture soil in wet and 
dry seasons 
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the background when only ambient air is used to 
fill the cell (5 cpm); a is the factor correcting gross 
counts to counts due only to decay of 222 Rn; CE 
is the efficiency of the cell-counter combination 
(cpm for a gas with known radon activity divided 
by the known activity; cell efficiencies ranged 
from 64% to 83%, averaging 76% efficient); 
V is the volume of the soil air sample; and the 
exponential term corrects for decay of 222 Rn 
(X = 1.258 x 10 ~ 4 min -1 ) in the cell between the 
time of collection (tfc) and the midpoint of the 
counting interval (tmp). Calculated errors for 
radon measurements reflect counting statistics and 
uncertainties in background and cell efficiency, 
and were typically about 10% of the measured 
activity. 

3.6. Laboratory measurements of radon production 

Radon production rates (reported as activity per 

gram of soil) were measured in the laboratory for 
both dry and wet soils. First, known quantities 
(50-100 g) of soil dried to constant weight were 
sealed in air-tight 100 ml jars. After waiting about 
14 days for the 222 Rn to come into secular equi- 
librium with the 226 Ra parent, air samples were 
removed by syringe and the activity of air in the 
jar was measured. After two or three measure- 
ments of the radon production rate from dry soil, 
distilled water was added to the jars to bring the 
gravimetric water content to typical field values 
(16-22% by weight), and measurements were 
repeated. 

3.7. Calculation of effective diffusivity in soil 

Total porosity (e) was calculated from measure- 
ments of bulk density (BD) and an assumed 
particle density (PD) of 2.65 g cm ~ 3 (Table 1 ): 

£=1-(BD/PD). (12) 

The mean particle density varies as a function 
of organic-C content and iron content, but the 
adjusted values would range only from 2.61 gem 1 
for the organic-rich surface soil to 2.73 gcm“ 3 at 
5 m depth where maximum total iron content is 
5.7% in these soils. Hence, these adjustments of 
particle density are trivial for calculations of 
porosity. 

Water-filled porosity (w) is simply the volu- 
metric soil water content determined from the 
calibrated TDR probes (see Subsection 3.2). 


Air-filled porosity is the difference between total 
; porosity and water-filled porosity (£ — w; Table 1 ). 
For calculation of effective diffusivity in soil (Z> 5 ) 
! using the nonaggregated model for dry media 
i (eq. (2)), only air-filled porosity (a) is needed. For 
| calculation of D s using the nonaggregated model 
! for wet media (eq. (4)), only total porosity (e) and 
» air-filled porosity (a) are needed. 

To use the aggregated soil model of eq. (7), the 
pore spaces must be partitioned into intra- and 
inter-aggregate spaces. Intra-aggregate porosity 
was estimated from volumetric soil water content 
at field capacity, which was measured by TDR 
probes about 48 h after a soaking rain near the 
end of the wet season in 1991. Field capacity is 
often estimated from water content measured at 
either —0.01 MPa or —0.03 MPa tension in the 
laboratory, but field capacity is more appro- 
; priately defined as the water content after a 
saturated soil has drained freely (Soil Science 
Society of America, 1987), which can result in a 
matric potential between —0.01 and —0.03 MPa, 
depending on soil texture (Papendick and 
Campbell, 1981). Our in situ estimates of field 
: capacity happened to agree reasonably well with 
laboratory estimates at —0.03 MPa tension, but 
no specific matric potential for field capacity is 
: assumed, and this approach conforms more closely 
to the formal definition of field capacity. It is 
assumed that, 48 h after rainfall has elevated the 
soil water content above field capacity, water has 
drained from the inter-aggregate pore spaces, 
leaving the intra-aggregate pore spaces 100% 
water-filled (£ ra = 6 fc ). Inter-aggregate porosity 
(£ er ) is estimated from the difference between total 
porosity and intra -aggregate porosity (c er = 
t — £ ra ; Table 1). For calculating diffusivity, it 
is assumed that the intra-aggregate spaces store 
water first and lose water last. When u'^c r# , £ cr is 
entirely air-filled so that a cr — c cr and a r& = £ ra — \\\ 
j When > £ ra , a ra = 0 and o er = £ — w. 

4. Results and discussion 

4. / . Radon profiles and estimates of diffusivity 

Activities of 222 Rn increased sharply from the 
surface to 1 m depth and then remained relatively 
constant from 1 m depth downward (Fig. 2). Depth 
profiles of 222 Rn such as these have been fitted 
t to exponential equations to estimate diffusivity 
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I 




222 Rn (kBq m 3 ) 


i 




222 Rn (kBq m' 3 ) 


Fig. 2. Depth profiles of radon activities in one forest soil and one pasture soil near the end of the wet season (May 
1992) and near the end of the dry season (November 1992). Circles show duplicate measured values, and error bars 
represent uncertainty due to analytical error (cell counting efficiencies, etc.). The solid lines show the predicted values 
using the nonaggregated model (eq. (4)), and the broken fines show the predicted values using the aggregated soil 
model (eq. (7)). 


(Dorr and Miinnich, 1990), but this approach 
assumes constant diffusi vity with depth. Signifi- 
cant gradients of soil water content in the upper 
two meters of the soil indicate that uniform 
diffusivity is unlikely, and calculated difTusi vi ties 
from both aggregated and non-aggregated models 
support this conclusion (Table 1 ). In addition, 
222 Rn activity (concentration) gradients in the 


upper part of the soil also depend on the radon 
activity at depth. Variations in the amount of 
air-filled versus water-filled pore space can cause 
significant changes in the observed deep soil radon 
activity on a seasonal timescale (Asher-Bolinder 
et a!., 1990; Washington and Rose, 1990). To 
account for these non-ideal conditions, we used a 
multi-box model that predicts 222 Rn activities 
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Fig. 3. Diagram of model used for predicting depth profiles of radon activities and C0 2 production, where a is the 
air phase, w is the water phase, C is concentration (or activity of radon) in the air phase, D is diffusivity, is the 
Bunsen coefficient, X is the decay constant for 222 Rn, and P is the radon production rate. The soil is divided into 50 
layers, each of 10 cm depth. 


throughout the soil profile using estimates of 
radon production and diffusion, and we compared 
these predictions with measured 222 Rn activities. 

To predict steady state radon activities in the 
soil atmosphere, we used the model illustrated 
in Fig. 3. The soil is divided into boxes (layers) 
of 10-cm depth, and the balance of radon is 
calculated through successive iterations using a 
10 minute time step until steady state is reached, 
usually in two weeks (4 half-lives of 222 Rn): 

d( TotRn,-)/d/ - P, - 2(TotRn,) 

-I- £>,_}( airRn,_ , — airRn,) 

— D,(airRn ( — airRn,+ j), (13) 

where for layer i, TotRn, is total 222 Rn (including 
air and water phases), airRn, is the 222 Rn activity 
in soil air, P t is 222 Rn production, D, is diffusivity. 


Production of 222 Rn was estimated in two ways: 
(i) the laboratory jar incubations described in 
Subsection 3.6; and (ii) by back -calculation of 
radon production from in situ radon activities in 
the deep (>3m) soil. In the laboratory jar 
method, dry soil radon production rates varied 
by depth, decreasing rapidly from values of 
13.5 Bq kg “ 1 (range: 9.5 to 1 7.0 Bq kg “ 1 ) in the 
top 10cm of soil to an average of 4.5 Bq kg -1 
(range: 2.7 to 8.8Bqkg -1 ) for depths of 50 to 
800 cm. The radon production rates were higher in 
wet soil, as has been observed by others (Stranden 
et al., 1984; Strong and Levins, 1982), with a mean 
of 19.3 Bqkg -1 for the top 10 cm, and a mean 
of 8.3Bqkg _1 (range: 4.5 to 11.7Bqkg “*) for 
depths greater than 50 cm. 

For the method based on calculations of 
measured in situ radon activities of deep soil, the 
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production rate (P)(Bq kg -1 soil) was calculated ; affected by diffusive loss of radon to the over- 
from the formula: lying air in addition to radioactive decay. The jar 

! incubation data clearly show that radon produc- 


A=Px BDx 


a + 0.22w* 


where A is radon activity (kBq m 3 air-filled pore . J 
space) for a given volume of soil; BD is the bulk \ 
density (g soil cm “ 3 bed spaced a the air-filled 
pore space (cm 3 air cm -3 bed space) and w is 
the water-filled pore space (cm 3 water cm" 3 bed 
space), and the Henry’s Law constant for 222 Rn 
is 0.22 (assuming instantaneous equilibration 
between air and water). This calculation is based 
on the assumption that the radon in the deep 
(>3m) soil is in secular equilibrium with the j 
226 Ra parent, which appears reasonable for soils 
below 2 m. Given a typical effective difTusivity of 
0.015 cm 2 s" 1 , the time required for radon to 
diffuse from 3 m to the surface is approximately 70 
days (square of the distance divided by the effective 
difTusivity), which is much longer than the half-life 
of 222 Rn (3.8 days). Hence, the equilibrium activity 
at depth is determined primarily by production 
and decay of 222 Rn, and not by loss to the surface. 
The only reason the deep soil radon content would 
not be at equilibrium would be if there were signi- 
ficant lateral exchange through the pit wall. The 
time required to diffuse the 1.7 m length of the 
sampling tubes would be about 25 days, again 
indicating that most of the radon would decay 
before escaping the soil via the pit wall. Analytical 
errors in the range of 10% of measured 222 Rn 
activities are greater than errors resulting from the 
influence of the pit wall. 

Radon production for deep soils calculated from 
eq. (14) averaged 12Bqkg"* in forests for both 
wet and dry seasons and averaged 16 Bq kg 1 in 
the pastures in the wet season and 12Bqkg -1 
in the dry season. These values are greater than 
production measured in either wet or dry jar 
incubations. We attribute this to the tendency of 
the wetted soils in the incubation jars to form 
clumps, which may have trapped radon within the 
soil. Thus, for most of the soil profile we chose to 
calculate the average radon production rates by 
the second method, which assumes equilibrium 
at depth. Radon production rates cannot be 
calculated from in situ measurements in the upper 
two meters of the soil, however, as these depths are 


tion rates in the top 50 cm of the soil are greater 
than those below, probably because radium is 
taken up by plants due to its similarity with other 
divalent cation nutrient elements and is con- 
centrated in the plants and plant detritus. We 
therefore used the wet jar incubation value 
(19 Bq kg -1 ) for the surface soil layer, and 
assumed a linear decrease to 50 cm, at which point 
the average deep in situ value (12 or 16 Bq kg" 1 , 
depending on the site and season) was assumed to 
be constant to 500 cm depth. Model results are not 
sensitive to changing the surface value by 60%. 

Seasonal variation in measured radon activities 
in forest soils (Fig. 2) can be explained using 
eq. (14) by the differences in water-filled porosity 
between wet and dry season and by assuming 
no change in production rate. In other words, 
a simple dilution effect caused by expansion and 
contraction of the air volume within the soil 
significantly affects the secular equilibrium value of 
radon activity when radon production is constant. 
Hence, variation in rates of radon production 
need not be invoked in this case to account for 
observed changes in radon activities. Larger 
seasonal changes in radon activities of pasture 
soils, however, require higher production rates 
to explain the high wet season deep soil radon 
activities. As we and others have observed in 
laboratory incubations (Stranden et al., 1984; 
Strong and Levins, 1982), radon production in 
soils is positively related to water content in soils, 
presumably due to reduced adsorption of radon 
onto surface particles when the surfaces are wet 
and because the soil water increases the emanation 
coefficient of radon. The 222 Rn generated by alpha- 
decay of 226 Ra possesses kinetic energy (alpha 
recoil energy), which is dissipated as the 222 Rn 
atom moves away from the site at which it was 
generated. In a fine grained medium, this distance 
may be great enough to transport the 222 Rn across 
air-filled pore spaces and into the next mineral 
grain. The presence of water in pore spaces serves 
to decrease the recoil range markedly, thus 
trapping the recoiled 222 Rn in pore spaces, and 
increasing the radon emanation observed in wet 
soils compared to dry soils (Nazaroff, 1992) The 
pasture soil w as very wet at depth in the wet season 
(Table 1), and perhaps these extreme conditions 
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Oil 


are needed for the water effect on radon production 
rates to be detectable by our in situ measurements. 

Returning to the model illustrated in Fig. 3, 
the soil properties shown in Table 1 were linearly 
interpolated between measured or calculated 
values to provide estimates for each 10 cm soil 
layer. Diffusivity was estimated for each soil layer 
using eq. (4) for the nonaggregated soil model and 
eq. (7) for aggregated soil model. Total radon in 
each 10-cm soil layer was distributed among air 
and dissolved phases using the Henry's law con- 
stant of 0.22 at 25 °C in fresh water (eq. (14)). We 
assumed instantaneous equilibration between soil 
air and water, that advective losses of radon were 
insignificant, and that diffusion gradients were 
based solely on activity (concentration) gradients 
in the soil air. The bottom soil box was assumed to 
have constant activity (equal to its secular equi- 
librium value). Predicted profiles of 222 Rn from the 
nonaggregated and aggregated models as well as 
observed 222 Rn activities are shown in Fig. 2. 

The choice between nonaggregated and aggre- 
gated soil models has no effect on modeled equi- 
librium activities of 222 Rn at depth (Fig. 2). As 
discussed above, the activity at depth is influenced 
primarily by rates of production and by the 
changing proportions of water-filled and air-filled 
pore space, rather than by diffusive losses to the 
surface or the pit wall. In contrast, the choice of 
diffusivity equations is important near the soil 
surface. The nonaggregated model significantly 
underestimates diffusivity in the top 1 m of soil, 
and this model predicts higher radon activities 
than were observed, including a large bulge in 
^Rn in the top 1 m of soil (Fig. 2). This failure of 
the nonaggregated model is not surprising, given 
that soils are aggregated media. Diffusivities 
calculated from the non-aggregated model are 
between 2 and 10 times lower that those from the 
aggregated model (Table 1 ). The aggregated soil 
model predicts 222 Rn activities that agree well with 
observed activities throughout the profile (Fig. 2). 
We concluded that the aggregated model ofeq. (7) 
(Millington and Shearer, 1971) works well for 
these soils and we adopted it for subsequent 
analyses of C0 2 production. 

4.2. Simplification of the aggregated media 
diffusivity model 

The disadvantages of the aggregated media 
model of Millington and Shearer (1971) are its 


seemingly complicated structure and its require- 
ment for differentiating among inter- and intra- 
aggregate porosity (eq. (7)). We have shown that 
the later can be overcome if data are available on 
volumetric water content at field capacity and if it 
is assumed that intra-aggregate spaces are water- 
filled and inter-aggregate spaces are air-filled at 
field capacity. 

Eqs. are not as difficult to work with 

as first impressions may indicate. Nevertheless, 
we have noted that the imposing nature of the 
aggregated media model could be mitigated by a 
simplification that is possible for the soils of this 
study and may be applicable elsewhere. The first 
term of eq. (7) accounts for diffusion through 
intra-aggregate pore spaces, but in all of the cases 
we applied it (all soil pits, all dates, and all depths), 
this term never accounted for more that 7 % of the 
effective diffusivity and usually only accounted for 
1 % or 2%. Hence, the first term could be dropped 
with little effect on the calculation of diffusivity, 
leaving the simplified formula: 

DJD 0 = a? T (aJe ct ) 2 . (15) 

It is likely that this simplified formula would 
perform satisfactorily in fine textured soils where 
water is often held tenaciously in small intra- 
aggregate pore spaces. The difference in volumetric 
water content, and hence in intra-aggregate air- 
filled pore space, between soils at field capacity 
and at —1.5 MPa matric potential is relatively 
small in fine textured soils. This simplified formula 
would be less likely to work in coarse textured 
soils, where the proportional change in water 
content is greater, but it deserves testing. 

Whether the full model or the simplified version 
is used, the aggregated media model is more 
appropriate for use in soils, especially in surface 
soils, than are other models. Gaseous diffusivity in 
soil can be estimated using only data on bulk 
density and water content at field capacity to 
characterize the soil and data on volumetric water 
content to account for the dynamics. 

4.3. CO 2 profiles and estimates of CO 2 production 
by depth 

Unlike 222 Rn, C0 2 does not decay, and the 
precision of our C0 2 concentration estimates is 
<2%, which is superior to the precision of the 
222 Rn measurements. Hence, the effect of mixing 
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Fig. 4. Horizontal profiles of C0 2 concentrations from 
the soil pit wall inward at 5 m depth in a forest soil 
measured early in the dry season (August 1993). The line 
shows the regression equation: 

^2<r) = 1 — e* - "*) + 2.5, 

where z is the tubing length (cm), C0 2(oo) is a fitted 
parameter for the asymptote C0 2 concentration (55.5 ml 
C0 2 air), a is a fitted exponential parameter, and 
2.5 ml C0 2 1 ~ 2 air is the measured C0 2 concentration in 
the air of the soil pit. 


of air from within the soil pit with air sampled 
near the end of the 180-cm long tubes is poten- 
tially within the range of detection of our C0 2 
measurements. Ideally, the tubes would be suf- 
ficiently long to avoid the pit wall effect for all 
gases, but this was not feasible because the longest 
dimension oTthe pit is only 180 cm, and an exten- 
sion had to^ taken on and off the auger for each 
auger bucket of soil removed beyond 180 cm. To 
test the effect of the pit wall on measured C0 2 con- 
centrations, we inserted tubes of 45, 90, 135, 180, 
and 270 cm length horizontally into the side of one 
of the forest soil pits at 5 m depth. An exponential 
fit of the horizontal C0 2 concentration gradient 
measured in August, 1993, early in the dry season, 
is shown in Fig. 4. The estimated asymptote 
concentration at infinite tubing length from an 
exponential nonlinear fit is 58 ml C0 2 I -1 soil air. 
The C0 2 concentration measured at 170 cm from 
the pit wall (about 10cm of the 180cm tube 
extends out into the pit) was 49 ml C0 2 1 1 soil 
air, or 85% of the concentration at the asymptote. 

This correction factor can also be independently 
estimated from estimates of diffusivity and using 
222 Rn as a proxy for C0 2 . We cannot use C0 2 
profiles directly because the equilibrium concen- 


tration of C0 2 is determined by production and 
diffusion, both of which are unknown. As already 
discussed, the 222 Rn concentration at equilibrium 
is determined by the rate of production and decay, 
and the decay rate of 222 Rn is known. The quotient 
of the 222 Rn activity at a given distance from the pit 
wall (^Rn.) divided by the 222 Rn activity at an 
infinite distance (^Rn^), where production and 
decay are at equilibrium, can be calculated from 
the exponential function described by Dorr and 
Munnich (1990): 


[ 222 RnL 

[ 222 Rn]oo 


=n-e-n 


(16) 


where z is the distance from the pit wall in cm, 
k — J{k/D 5 ), k is the decay constant for 222 Rn 
(2.0978 x 10“ 6 s~*), and D s is the effective dif- 
fusivity of 222 Rn in soil calculated from eq.(7). 
Constant diffusivity and constant 222 Rn produc- 
tion along the horizontal length of the sampling 
tube are assumed. This quotient must then be mul- 
tiplied by 1.2 to account for more rapid diffusion of 
C0 2 relative to 222 Rn. The quotient from eq. (16), 
adjusted for C0 2 , is 0.85 for 170 cm distance from 
i the pit wall at 5 m depth and measured soil water 
contents of August, 1993, and using the diffusivity 
estimates from the aggregated media model 
(eq. (7)). The excellent agreement between the 
: quotient calculated from diffusivity estimates and 
the correction factor estimated from the measured 
horizontal C0 2 profile shown in Fig. 4 indicates 
that the aggregated media model yields reasonable 
estimates and that the measured CO : concentra- 
tions within the deep soil can be reliably adjusted 
to account for diffusional losses through the pit 
walls. The quotients calculated from eq. (16) 
ranged from 0.80 to 0.95 for the data from our 6 
soil pits in pastures and forests and during wet and 
dry seasons. 

Depth profiles of measured concentrations of 
C0 2 (not corrected for the pit wall effect) are 
shown in Fig. 5, and they increased with depth to 
at least 8 m. Concentrations at depth are higher in 
forest soils than in pasture soils. Concentration 
alone, however, does not indicate the strength of 
C0 2 production sources, because high concentra- 
tions can result from either high rates of produc- 
tion or low rates of diffusion. 

To estimate C0 2 production, we first solved 
eq. (16) for the asymptote C0 2 concentration at 
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Fig. 5. Depth profiles of C0 2 concentration (means of duplicate analyses) in forest and pasture during the wet season 
(squares) and the dry season (circles). Data are from the same sites and dates as shown for radon activities in Fig. 2. 


infinite tube length, using estimates of diffusivity 
from the aggregated media model (eq. (7)) and 
measured concentrations at 170 cm tube length at 
depths below 1 m. We then applied the corrected 
C0 2 concentration profiles and the diffusivity 
estimates to a model based on Fick’s first law. The 
C0 2 production (Pco 2w ) within each 10-cm soil 
layer of the model represented by Fig. 3 was 
estimated from the difference between the C0 2 flux 
calculated from Fick’s first law at the top and 
bottom of each layer; 

Pc o 2 .,= [Ax(C,- C,_,)/z] 

-t A + ,x(C (+I -C,)/z] (17) 

where for layer i, Z), is the diffusivity calculated 
for the aggregated medium model (eq. (7)), C, is 
the soil C0 2 concentration, z and is the depth 
increment (10cm). Concentrations of C0 2 were 
linearly interpolated to estimate C0 2 concentra- 
tions within each 10 cm layer. Because of uncer- 
tainties in interpolations of C0 2 concentrations 
and soil water contents, the C0 2 production 
estimates for individual 10 cm layers may not be 
reliable, but reliability increases as estimates from 
the 10 cm layers are summed for larger depth 
increments that include measured values of input 
variables. We also suspect that this approach 
underestimates C0 2 production in the layers 
above 1 m, because C0 2 concentrations are not 
linear or unimodal in this region, but rather 
exhibit “hot spots” of high C0 2 production and 


concentration. Finally, C0 2 production in the 
forest floor litter cannot be estimated directly by 
eq. (17) because we do not have estimates of 
diffusivity or measurements of C0 2 concentration 
in the forest floor layer. 

The ecological question of interest is the relative 
importance of surface and deep soils as sources of 
soil CO 2 in forests and pastures. For this purpose, 
we have grouped the production estimates for 
three large depth intervals (Table 2). Production 
below 5 m is simply the flux calculated from Fick’s 
first law at 5 m depth — any C0 2 diffusing upward 
at that depth is assumed to have been produced 
below. Production below 5 m depth may be under- 
estimated if there is also a downward flux of C0 2 
at some depth, either through the soil atmosphere 
or dissolved in drainage water. We cannot address 
at this time the magnitude of export of C0 2 via the 
groundwater. The production estimate for the 1 to 
5 m depth increment is the sum of the C0 2 produc- 
tion rates calculated from eq. (17) for all of the 
10 cm soil layers between 1 and 5 m. The estimate 
for the litter layer plus the top 1 m of mineral soil 
is the difference between measured fluxes of 
C0 2 from the soil/litter surface (Subsection 3.3, 
Table 2), and the Fick’s law flux calculated at 1 m. 
These estimates were calculated for each of the 
three soil pits in forest and pasture, and the means 
are reported in Table 2. 

Surface fluxes of C0 2 from the forest are about 
twice the rates measured in the degraded pastures 
(Tabic 2). The forests are more productive eco- 
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Table 2 Summary of measured CO, fluxes from the soil surface and calculated rates of production of CO 2 
™ k i: f son Ji litter (L): means (rounded to the nearest JO mg) and standard errors 


are of three soil pits per ecosystem type 


Wet season (May 1992) 


Dry season (November 1992) 


forest 


pasture 


forest 

pasture 


mean 

SD 

mean 

SD 

mean 

SD 

mean SD 




mg C0 2 -C m 

- 2 h-' 



surface flux 
production: L to 1 m 
production: 1 to 5 m 
production: below 5 m 

290 

240 

30 

20 

20 

24 

5 

7 

140 

100 

30 

<10 

6 

16 

11 

2 

240 

200 

10 

30 

20 

26 

15 

15 

110 23 

70 26 

20 4 

10 0 


systems than are degraded pastures, and below- 
ground inputs of C are apparently higher in the 
forests than in the pastures. About 70-80% of the 
production in both ecosystems occurs within the 
top 1 m of soil, including the forest litter layer. 
Very little C0 2 ( *5 10 mg C m ' 2 hr ' 1 ) is produced 
below 5 m depth in the pasture. Production of C0 2 
below 5 m depth is also low in the forest sites 
relative to the surface, but the deep soil C0 2 
production rates in forests are more than double 
those in the pastures. 

4.4. Ecological significance of research findings 
The deep roots of eastern Amazonian forests 
enable the vegetation to extract water stored deep 
in the soil and thus retain green foliage during the 
long dry season (Nepstad et al., 1994). In contrast, 
pasture vegetation turns brown during the dry 
season and pasture vegetation extracts less water 
from deep soils. Active extraction of soil water by 
deep forest roots must be accompanied by root 
respiration and microbial decomposition of root 
inputs and by concomitant production of C0 2 . 
Although most of the C0 2 production is concen- 
trated in the top 1 m of soil, as are most of the 
roots (Nepstad et al., 1994), the deep soil processes 
are not trivial. Our estimates of C0 2 production 
by depth (Table 2) are consistent with the depth 
distribution of fine root biomass reported by 
Nepstad et al. (1994), in which 70% to 80% 
occurred above 1 m and the remainder below 1 m 
depth. Changes in land use, from deep-rooted 
forest vegetation to more shallowly rooted pasture 
vegetation, could change the depth distribution of 


root biomass, C0 2 production, and the dynamics 
of soil organic matter, ir half of the C0 2 produced 
below 1 m depth emanates from microbial decom- 
position (the other half emanating from root 
respiration), the annual rate of microbial decom- 
position of soil organic matter below 1 m depth 
would be about 0.2 kg Cm' 2 , which indicates 
. that, on average, about 1 % of 20 kg C m * stock 
of soil C found between 1 m and 8 m depths turns 
over annually. If C inputs to deep soils decline 
when forests are converted to pastures, the deep 
soil C stocks could change significantly over 
decades to centuries. Fully understanding the role 
of land use change on the global carbon cycle will 
require consideration of these deep soil processes. 
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